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We provide a toolbox for continuous variables quantum state engineering and characterization 
of biphoton states produced by spontaneous parametric down conversion in a transverse pump 
conhguration. We show that the control of the pump beam’s incidence spot and angle corresponds to 
phase space displacements of conjugate collective continuous variables of the biphoton. In particular, 
we illustrate with numerical simulations on a semiconductor device how this technique can be used 
to engineer and characterize arbitrary states of the frequency and time degrees of freedom. 
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Spontaneous parametric down conversion (SPDC) ex¬ 
periments play a prominent role in the field of quantum 
information and communications. Single and multiple 
photon pairs generated through SPDC display entangle¬ 
ment in multiple degrees of freedom (DOF) which are 
fundamentally different from each other. When com¬ 
bined or independently accessed, they constitute a power¬ 
ful platform for experimental demonstrations of quantum 
protocols. Discrete DOF, such as polarization and orbital 
angular momentum, are currently used to implement 
quantum logic gates and protocols [1], test the non-local 
properties of quantum mechanics [2-5] , and realize quan¬ 
tum cryptography [6-8] and teleportation [9, 10]. DOF 
associated to observables with a continuous spectrum, 
like the quadratures of the electromagnetic field, poten¬ 
tially offer the same versatility as discrete ones in the field 
of quantum information. Continuous DOF in the single 
photon regime, such as frequency, transverse momentum 
or position display a perfect analogy with a multi-photon 
single mode continuous variables (CV) [11, 12]. Conse¬ 
quently, they constitute an attractive platform to real¬ 
ize CV quantum information protocols [13] that are usu¬ 
ally associated to the single mode multi-photon config¬ 
uration. For example, an appealing aspect of using sin¬ 
gle photon’s transverse coordinates in this field is their 
relatively easy manipulation with readily available opti¬ 
cal devices, such as spatial light modulators (SLMs) and 
lenses [12], circumventing the difficulties encountered in 
multi-photon CV strategies to implement non-gaussian 
operations, which are essential ingredients of universal 
quantum computation [13]. For these reasons, the study 
of entanglement in CV in the single photon regime is 
a valuable strategy to demonstrate CV-based quantum 
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gates [14], quantum key distribution [15, 16], error cor¬ 
recting codes [17, 18], quantum metrology protocols [19] 
and to study quantum to classical transitions [20, 21]. 
Finally, we can mention that continuous DOF can be 
combined to discrete ones to implement conditional op¬ 
erations [22-24]. 

CV entanglement in photon pairs can be generated 
in different DOF via SPDC; one example are the spa¬ 
tial transverse DOF of photon pairs produced in non¬ 
linear bulk crystals. Tayloring the spatial properties 
of the pump beam has been previously employed to 
engineer [25, 26] and detect quantum states of photon 
pairs with different properties : in [27] it was shown 
that controlling the transverse phase properties of the 
pump beam could lead to bunching or anti-bunching of 
a photon pair. Refs. [22, 28] demonstrate that the use 
of SLMs put in the arms of an interferometer can lead 
to the measurement of arbitrary momenta of the quan¬ 
tum state associated to the transverse photonic degrees 
of freedom. 

Using frequency as a CV degree of freedom instead 
of the transverse coordinates presents a series of poten¬ 
tial advantages: in general, the transverse variables of a 
propagating field are entangled in the two-dimensional 
orthogonal spatial coordinates x and y, or equivalently, 
the momentum coordinates Px and Py] this renders arbi¬ 
trary state production and measurement more challeng¬ 
ing, and limits the applications of such degrees of freedom 
in quantum information tasks, since the quantum state 
associated to each spatial direction is entangled and can 
be cross-correlated between x and y directions. More¬ 
over, this spatial entanglement is jeopardized when cou¬ 
pling into optical fibers or waveguides, while frequency 
states remains robust in such devices. For all these rea¬ 
sons, we will illustrate our ideas focusing on the spectral 
degree of freedom. Indeed, the characterization of arbi¬ 
trary frequency entangled states of photon pairs draws 
a lot of attention as illustrated by recent works on fre¬ 
quency dependent intensity correlation measurements of 
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the photon pair [29] , the reconstruction of biphoton wave- 
functions generated with a monochromatic pump [30] and 
amplitude sensitive tomography techniques in the time- 
energy space [31]. 

In the waveguided regime, frequency correlated, un¬ 
correlated and anti-correlated photon pairs can be pro¬ 
duced by modifying the spectral properties of the pump 
beam [32, 33]. Nevertheless, even if the nature of corre¬ 
lations can be controlled, pump bandwidth modification 
alone cannot lead to arbitrary CV state engineering, an 
essential tool for CV quantum information and commu¬ 
nications, fundamental tests of quantum mechanics and 
quantum metrology. 

In this manuscript, we provide a toolbox to exploit 
the capabilities of SPDC to create, measure and char¬ 
acterize entanglement in the CV DOF of a photon pair. 
Our strategy is based on the combination of pump spa¬ 
tial engineering with the possibility of characterizing a 
Wigner function of the photon pair using a Hong-Ou- 
Mandel (HOM) [34] type experiment [35]. 

In order to clarify these concepts, we present the case 
study of a tranversally pumped semiconductor waveguide 
[36, 37] potentially having a great versatility in the con¬ 
trol of the biphoton frequency correlations [29]. Indeed 
it has been shown that the pumping configuration allows 
one to modify the phase matching conditions by simply 
changing the spatial properties of the pump beam [38-40], 
e.g. its waist or angle of incidence. However, more com¬ 
plex quantum state engineering has not been explored so 
far. 

The properties of the photon pair are described by the 
biphoton wavefunction which takes a complex form due 
to the constraints imposed by momentum and energy 
conservation. By neglecting group velocity dispersion 
(which is justified in the spectral range of the generated 
photon pairs) and in the narrow bandwidth limit for the 
pump beam, we can write the state of the pair under the 
following form: 

|^')=Xr JJ dLOsduJif+{uJs+uJi)f-{cos-uJi)\cos,coi) (1) 

where represent the signal and idler frequencies, re¬ 
spectively, and xr is a normalization constant [41]. The 
function /+ corresponds to the spectrum of the pump 
beam reflecting the condition of energy conservation 
{ujp = ujs uJi = The phase of this function is 

affected by the presence of a microcavity formed by two 
Bragg mirrors as can be seen in the non flat phase profile 
in Fig. Ic. The narrow band assumption allows us to con¬ 
sider the phase-matching dependent part of the biphoton 
state as a strictly antidiagonal function /_: 

r^/2 ^ (ujs-uji)z 

f-{uj- = uJs - = / dz(p{z)e ^9 , (2) 

J-LI2 

where L is the length of the sample, Vg is the average 
group velocity of the signal and idler photons at fre¬ 
quency degeneracy and ip{z) = The func¬ 

tion ^{z) is the spatial profile of the pump beam and 
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-, where Ug i corre- 
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spends to the signal and idler effective refractive indices 
at degeneracy, ujp^^ is the central frequency of the pump 
beam and 6>deg is the pumping incidence angle for which 
degeneracy occurs. 

If the dimensions of the waveguide are large with re¬ 
spect to the pump waist, i.e. the limit where L ^ oc, 
/_ can be approximated by the Fourier Transform of the 
spatial profile of the pump beam: 


We start by considering the situation depicted in 
Fig. la where a gaussian pump beam with waist Wp 
impinges onto the source at an angle 0 and position 
zq. The field distribution along the axis ^ is ^{z) oc 
^-{z-zo) cos eiwp^i(ksin 9)z therefore /_ reads: 


/_(cc;_) oc (4) 

with To = z^jvg^ Auj = Vgl2wp and = {ksinO — 
^deg)'^^ ~ SOvgUOpjc [42]. Fig. lb and Fig. Ic represent 
the biphoton wavefunction norm and phase distribution 
numerically simulated for a pump impinging at the de¬ 
generacy angle 6>deg at zq = 0 with waist Wp = 200 pm ^ 
L = 2mm, central wavelength Xp = 775nm and pulse 
duration Tp = 3.2ps [43]. 

Instead of the complex-valued biphoton wavefunction, 
a more convenient representation is given by chronocyclic 
Wigner functions which are the time-energy analogs of 
the phase-space Wigner functions. With this approach, 
time and energy properties of the single photons of the 
pairs are illustrated with real quantities [44, 45], the same 
way pulses of light are depicted in the domain of ultra¬ 
fast optics [46]. In this work, we describe, instead of the 
properties of the isolated photons of the pair, their cor¬ 
relations along the antidiagonal part /_ of the biphoton 
wavefunction. This represents the quasi-probability dis¬ 
tribution of the biphoton as a function of the detuning Q 
between signal and idler frequencies and the time delay 
conjugated variable r. The corresponding Wigner func¬ 
tion W- is given by: 

/ c» 

dco’_/_(r^ — + . (5) 

-oo 

Using the expression obtained in (4), we see that this cor¬ 
responds to a gaussian Wigner function centered at point 
r = To, U = of widths AQ = 1/Ar = Auj = Vgj^Wp 
that is equivalent to the representation of a coherent state 
(see Fig. Id). Thus, in this situation, shifting the pump¬ 
ing spot zq is equivalent to realizing displacements in the 
r axis of the phase space while changing the angle of in¬ 
cidence 0 of the pump beam corresponds to shifting the 
state along the fl axis. 
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FIG. 1: (Color online) (a) Counter-propagating phase-matching scheme implemented in a semiconductor 
microcavity waveguide. The characteristics of the pump pulses allow to tune the time-energy properties of the 
biphoton, (b) Norm and (c) phase of the biphoton wavefunction for Wp = 200p,m, Xp = 775nm, Tp = 3.2ps and 
0 = 6>deg- The phase pattern in (c) is due to the resonance of the pump in the microcavity, (d) Corresponding 

Wigner function. The ft axis is given in units A = 



FIG. 2: (Color online) (a) Pump illumination scheme to generate a compass state : two pairs of beams impinge onto 
the waveguide at two spots Za and equidistant from the center of the source, each pair consisting of two beams 
tilted symmetrically with respect to the degeneracy angle {0a,b = ^deg i (b) Norm and (c) phase of the biphoton 
wavefunction of the corresponding cat state with 60 = 9.37' ^ \za — Zb\ = 1mm and for each beam Wp = 200pm, 

Xp = 775nm, Tp = 3.2ps and 0 = ^deg- (b) Corresponding Wigner function with AA = 1.37nm and Ar = 10.8ps. The 

ft axis is given in units A = ^^ft. 


More complex states can be obtained by engineering 
the pump beam. Indeed two identical beams imping¬ 
ing at Za and z^ will generate a superposition of two co¬ 
herent states displaced along the r axis. In the limit 
\za — Zb\ ^ such states are almost orthogonal, rep¬ 
resenting a superposition of two distinct quasi-classical 
states (Schrodinger cat-like states). An analog superpo¬ 
sition is obtained along axis ft by using 2 different angles 
of incidence Oa and Ob, impinging at the same point zq. 
Quasi-orthogonality is obtained for \0a — 0b\ ^ —. 

We can generalize these Schrodinger cat-like states us¬ 


ing more complex configurations of pump beams. As an 
illustration, we choose a compass state (see Fig. 2), a 
superposition of four coherent states presenting interest¬ 
ing applications in quantum metrology, as pointed out in 
[47, 48]. In order to obtain such state, a set of 4 different 
pump beam configurations is required: 2 pairs of beams 
impinging at 2 different points separated by a distance 
A 2 :, each pair consisting in 2 beams tilted symmetrically 
with respect to the degeneracy angle as shown in Fig. 2a. 

The procedure detailed above to generate quantum 
states with different properties and applications can be 
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FIG. 3: Two possible strategies to realize the modified 
HOM experiment leading to the measurement of the 
Wigner function. Left: displacements in time and 
frequency are realized in each arm of the interferometer 
after the production of the photon pair. Right: in a 
completely analogous set-up, displacements in time and 
frequency are realized by pump engineering. While 
modifying the incidence angle of the pump beam, and 
consequently, the phase matching condition, 
corresponds to displacements in frequency, modifying 
the pump’s incidence spot corresponds to displacements 
in time. 


generalized, leading to the creation of arbitrary continu¬ 
ous variables quantum states, since any state can be con¬ 
structed from the over-complete basis of coherent states. 

We now discuss how pump engineering can be used to 
characterize arbitrary CV states using the direct mea¬ 
surement of the Wigner function in all points of phase 
space. In [35], it was shown how a generalization of the 
HOM experiment leads to the direct measurement the 
Wigner function in the ^4, r phase space. Such gen¬ 
eralization consists of considering displacements in the 
frequency degree of freedom of the photons, in addition 
to time displacements. Of course, displacing either one 
of these parameters modifies the distinguishability be¬ 
tween the photons in each arm of the HOM interferome¬ 
ter, as depicted in Figure 3. Time displacements can be 
realized relatively straightforwardly by simply modifying 
the optical path in each arm of the HOM interferometer. 
However, the frequency displacements required to recon¬ 
struct the Wigner function are quite broad with respect 
to the performances of currently available optical modu¬ 
lators [49]. Pump engineering provides an alternative so¬ 
lution to realize both time and frequency displacements, 
dramatically simplifying the direct measurement of the 
Wigner function and the CV state characterization. 

As discussed previously, modification of the pump 
beam’s incidence angle corresponds to displacing the cen¬ 
tral frequency of the symmetric part of the wave function 


associated to the photon pairs, while modification on the 
incidence point corresponds to time displacements. Us¬ 
ing these ingredients, one can devise a procedure for the 
complete Wigner function measurement as follows: in a 
first step, an initial state is engineered. Such state is the 
one to be characterized. Running the HOM experiment 
with no frequency or time displacement leads directly to 
the value of IU_(0,0), the Wigner function at the ori¬ 
gin [35]. Then, tilting the incident pump beams by a 
given amount, and repeating the same HOM experiment 
is equivalent to displacing in frequency the original state 
and measuring its Wigner function, leading to the value 
of IU_(U,0). Analogously, displacing the pump beam in 
the 2 : axis and repeating the HOM strategy leads to the 
value of IU_(0,r). It is clear that, by combining differ¬ 
ent tilting angles and displacements, one can obtain an 
arbitrary point of the Wigner function and reconstruct 
IU_(U, r) for all the values of U and r. In order to char¬ 
acterize the quantum state, the magnitude of the dis¬ 
placements in both axis of phase space should cover the 
region where |/_(co’_)p has a significant value. This cor¬ 
responds to realizing the pump’s angular displacements 
in an interval ~ ^6 = 18.7', while its imping¬ 
ing position is displaced of /\z ^ Za — = 1mm (see 

Fig. 2). Notice that, as shown in [50] the proposed strat¬ 
egy presents the advantage of not being limited by the 
detector’s response time for measuring highly oscillating 
fringes or phase space structures associated with sub- 
Planck scales [47]. Note that the resolution required for 
displacements along the 2 :-axis is easily achieved; as far 
as the angular displacement is concerned, since AO is of 
the order of 2.7 mrad, a resolution of 100 prad is suffi¬ 
cient to resolve the fringes. This is also achievable with 
thermal stabilization and stable mechanical mountings. 

One may argue that modifying the pump, in reality, 
modifies the state to be measured instead of the measur¬ 
ing apparatus that is probing different points of the phase 
space; but this is common practice in quantum measure¬ 
ment strategies, where the modification of the settings 
of the measurement apparatus is formally equivalent to 
that of the state-to-be-measured. This equivalency has 
been previously used, for instance in the context of cav¬ 
ity quantum electrodynamics in [21]. There, the Wigner 
function is directly measured using a Rydberg atom inter¬ 
acting dispersively for a fixed time interval with the field 
of a high quality micro wave cavity. The set-up is kept 
the same for all points of phase space and the quantum 
state of the field in the cavity to be measured is displaced, 
and consequently modified, through the application of a 
coherent field. Also in [30], the use of a monochromatic 
pump generating pairs of photons through SPDC is com¬ 
bined with modification of the pump frequency (or shifts 
in the temperature of the crystal) to circumvent the dif¬ 
ficulty of broad frequency displacements. In the present 
case, changing the pump spatial configuration is equiva¬ 
lent to displacing the state to be measured. 

In conclusion, we have shown how spatial pump en¬ 
gineering can be used to generate arbitrary CV fre- 
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quency states of a photon pair and directly characterize 
it through a HOM-like experiment. The combination of 
the variety of pumping geometries with that of the pos¬ 
sible design of integrated photonic circuits [51] will allow 
the realization of complex and versatile quantum pho¬ 


tonic chips [40]. The generalization of this technique to 
prepare quantum states with higher photon numbers, as 
well as to simulate the dynamics of CV states and the 
application of different quantum operations to it is an 
interesting open perspective that will be investigated. 
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